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ABSTRACT 

In this paper we investigate the X-ray spectra and X-ray spectral variabihty of compact X-ray 
sources for 3 Chandra observations of the Local Group galaxy M33. The observations are centered on 
the nucleus and the star forming region NGC 604. In the observations 261 sources have been detected. 
For a total of 43 sources the number of net counts is above 100, sufficient for a more detailed spectral 
fitting. Of these sources, 25 have been observed in more than one observation, allowing the study 
of spectral variability on ^months timescales. A quarter of the sources are found to be variable 
between observations. However, except for two foreground sources, no source is variable within any 
observation above the 99% confidence level. Only six sources show significant spectral variability 
between observations. A comparison of Nh values with HI observations shows that X-ray absorption 
values are consistent with Galactic X-ray binaries and most sources in M33 are intrinsically absorbed. 
The pattern of variability and the spectral parameters of these sources are consistent with the M33 
X-ray source population being dominated by X-ray binaries: Two thirds of the 43 bright sources 
have spectral and timing properties consistent with X-ray binaries; we also find two candidates for 
super-soft sources and two candidates for quasi-soft sources. 

Subject headings: X-rays: binaries — galaxies: Local Group — galaxies: individual (M33) 



1. INTRODUCTION 

M33 is a late-type spiral galaxy, Sc II-III, and the third 
largest galaxy in the Local Group. It is a unique galaxy 
in the Local Group since morphologically it is of inter- 
mediate type between the large early-type spiral galaxies 
and the numerous dwarf galaxies. Other galaxies of this 
type cannot be investigated with the same depth even 
with Chandra . At a distance of 840 kpc (Freedman 
et al. 1991) from the Milky Way (MW) M33 is the sec- 
ond nearest major galaxy. It spans roughly 73'x45' on 
the sky. The line-of-sight absorption column density is 
small, TVh - 6 X 10^" cm^^ (Stark et al. 1992). M33 is 
more actively star forming than either the MW or M31 
(Hippelein et al. 2003), particularly compared to its much 
smaller mass. 

M33 has been studied with every X-ray mission since 
Einstein (Markert & Rallis 1983). But only recently have 
high angular resolution and high sensitivity instruments 
like Chandra and XMM allowed us to study the X-ray 
source population in depth. Grimm et al. (2005) and 
Pietsch et al. (2005) have provided source lists and fluxes 
for the X-ray source population in M33, from the Chan- 
dra and XMM observations respectively. In this paper 
we follow up on the Grimm et al. (2005) Chandra survey 
and present an analysis of the X-ray spectra and variabil- 
ity behavior of X-ray sources in M33. Apart from fluxes, 
variability and spectral energy distributions are an im- 
portant diagnostic tool for understanding the emission 
mechanism(s) in X-ray sources, and for classifying these 
sources. Moreover, spectral and time variability are the 
main characteristics of X-ray binaries, as shown by the 
detailed work done for X-ray binaries in the Milky Way, 
see e.g. van der Khs (2005). 

The paper is organized as follows. In Sec. 2 we describe 
the data processing, followed in Sec. 3 by the analysis 
procedures. The results and their implications for the na- 
ture of the X-ray source population in M33 are discussed 



in Sec. 4; we also discuss individual bright sources in 
detail in Sec. 4.4. We conclude with a summary of the 
results in Sec. 5. The figures showing the results for the 
spectral analysis for sources with more than one observa- 
tion are shown in the Appendix in the electronic edition 
only 

2. DATA 

M33 has been observed with the ACIS instrument on 
Chandra four times, see Grimm et al. (2005). In this pa- 
per we use three observations whose Obslds and dates 
are given in Table 1. Due to its angular extent the 
observed parts of M33 cover all active chips, the stan- 
dard ACIS-S configuration for Obsid 786, and the stan- 
dard ACIS-I configuration for ObsId 1730 and ObsId 
2023. There is considerable overlap between the differ- 
ent observations. However, due to the decreasing reso- 
lution/sensitivity with increasing off- axis angle only the 
inner part of M33 (^ 8-10 arc minutes) has a significant 
number of sources in two observations, 786 and 1730. A 
fourth observation, ObsId 787, was disregarded because 
it was aimed at studying the nucleus and suffers from 
both high background and small FOV. 



TABLE 1 
List of ACIS observations of M33 



ObsId 


Date 


Aim point 


Duration [ks] 


786 


2000-08-30 


Nucleus 


45 


1730 


2000-07-12 


Nucleus 


45 


2023 


2001-07-06 


NGC 604 


90 



The data from the three observations were processed 
according to the standard data processing procedure 
with CIAO versions 3.1, including exposure correction. 
Source detection was performed with wavdetect with 
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scales of 1, 2, 4, 8, 10, 12, and 16 in tlie energy range 0.3 
8 keV. Tire signal detection threshold was set to 10^®. 
Source regions correspond to the 95% encircled energy 
area. For more details of the data analysis see Grimm 
et al. (2005). 

3. ANALYSIS 

We separate the analysis in three parts, short term, 
long term, and spectral variability. The different meth- 
ods are discussed in the following subsections. We as- 
sume that the variability of individual sources is inde- 
pendent, so there is no correlation between short term 
and long term variability. Because the number of counts 
observed in a single observation is not very large even for 
bright sources we restrict the spectral variability analy- 
sis to a comparison between different observations, when 
applicable. 

The nucleus was excluded from the following analysis 
because it suffers from strong pile-up in two observations. 

3.1. Short and long term variability 

In order to establish short timescale variability we per- 
form a Bayesian block analysis of the lightcurves of indi- 
vidual sources. A Bayesian block analysis computes the 
best approximation to the lightcurve shape in terms of 
piccewisc constant flux levels or blocks. The discrimina- 
tor between a single flux level for the whole lightcurve 
or two flux levels is the ratio of likelihoods of describing 
a data segment with one or two blocks. In this analysis 
the algorithm used is iterative. It starts with the whole 
lightcurve and subsequently divides it until the likelihood 
ratio for dividing a lightcurve segment becomes smaller 
than a predefined prior, in this case corresponding to 
-^99% confidence level. The Bayesian block analysis is 
particularly suited for burst-like variability. For more 
details about the principle of Bayesian blocks and the 
algorithm used in this analysis see Scargle (1998). The 
implementation used here is the same as used by the 
CHAMP project and described in (Kim et al. 2004). 

Since a Bayesian block analysis uses only the photon 
arrival times for computation of the blocks, no binning 
is necessary. Therefore there is no intrinsic restriction 
to the number of photons the algorithm is applicable to. 
Obviously establishing variability with a certain confi- 
dence becomes less likely for fainter sources. 

Because a Bayesian block analysis is particularly sen- 
sitive to burst-like variability we also performed a search 
for periodic variability with the XRONOS v5.21 tool ef- 
search. This analysis did not yield any source with sig- 
nificant short term variability. 

To establish long timescale variability we use simple 
Poisson statistics. This is justified because, with the 
exception of two foreground sources, no source exhibits 
strong short term variability. We compare the difi^erence 
in fluxes between two observations to the quadratically 
added errors for the fluxes. If a source is undetected in an 
observation in which it was in the fleld of view of Chan- 
dra , we calculate; an upper limit to the source flux using 
an algorithm developed by Kraft et al. (1991). Since the 
upper limit value already corresponds to the 99% con- 
fidence level, for a comparison with a detected flux we 
compare the 3 sigma error of the detection directly with 
the upper limit. 



3.2. Spectra and spectral variability 

For 43 out of the total 261 sources the number of counts 
in at least one observation was larger than 100 net counts, 
which wc consider sufficient to attempt spectral fitting. 
Of these 43 sources 25 have been observed in at least 
two observations. For these sources we also compare the 
spectral properties with time. The spectral fitting was 
done with XSPEC vll.3.1. Because of the generally low 
number of counts the fitting was done for all sources with 
Cash statistics (Cash 1979). At 100 counts we expect 
only ~3 background sources based on the CDF-N LogN- 
LogS (Alexander et al. 2003), so the vast majority of 
these sources are likely to belong to M33. 

We first fit the bright sources with two simple absorbed 
spectral models, power law and bremsstrahlung, to check 
the validity of our assumptions about the spectral shape 
of the faint sources that do not have sufficient counts 
for detailed spectral modeling (see derivation of the X- 
ray luminosity function in Grimm et al. (2005)). The 
absorption was in one case a free parameter in the other 
case it was fixed to the Galactic value. After validating 
that the assumptions about the general power law shape 
are correct, we proceeded to fit more complex model to 
the bright sources. 

The results of the simple spectral fitting are pre- 
sented in Sec. 4. Spectra of sources with more than 
one observation with best fit values and correspond- 
ing confidence contour plots are shown in Appendix 
B in the electronic edition, except for the spectra of 
M33 X-7 (CXO J013334.1+303210) and M33 X-9 (CXO 
J013358.8+305004), which are shown in Figs. 7 and 8, 
and will be discussed in Sec. 4.4. 

Based on the confidence regions for the fit parame- 
ters, shown in Appendix B, only one of the X-ray sources 
shows significant spectral variability from observation to 
observation, M33 X-4 (CXO J013315. 1-^305317). For the 
other sources in different observations none of the param- 
eters of a source spectrum show changes corresponding 
to more than 99% confidence level. This is despite the 
fact that some sources show significant time variability 
between observations. However, the number of counts 
available for most sources for fitting is insufficient to de- 
termine spectral parameters to an accuracy good enough 
to compare two observations. Moreover, degeneracies be- 
tween model parameters complicate the establishing of 
variability. 

Another approach to spectral variability is the use of 
hardness ratios. These are relatively crude estimators 
of spectral changes but because of the smaller number 
of degrees of freedom compared to a spectral fit can be 
statistically preferable. Wc therefore take all sources 
that were observed in at least two observations, divide 
each observation in three equally long parts, and con- 
struct hardness ratios from these intervals. Because only 
sources with more than 100 counts ar(^ usc^d in this part 
of the analysis the separation in three parts still gives suf- 
ficient number of counts (> 20 — 30) in each time bin for 
a hardness ratio analysis (Prestwich et al. 2003). Note 
that Obsid 2023 is roughly twice as long as the other 
observations. The energy bands used are the same as in 
Grimm et al. (2005); 0.3-1.0 keV for the soft band, 1.0- 
2.1 keV for the medium band, and 2.1-8.0 keV for the 
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hard band. The hardness ratios are defined as 

where 5, M, and H are the background and exposure 
corrected counts in the soft, medium, and hard band, 
and T arc the corrected counts in the whole energy band. 
To compare two observations it is important to take into 
account the different effective areas of the source location 
in each observation, especially if the source is located on 
a front-illuminated chip and in another observation on a 
back-illuminated chip. We compute the effective area of 
the source in each observation for each energy band. For 
comparison with the other observations of the source we 
normalize the effective area in each band by the effective 
area value of the aim point of Obsid 1730. Note that the 
choice of the normalization constant is arbitrary. The 
results arc combined for all available observations. The 
results for all sources that show evidence of variability 
are shown in Appendix B in the electronic edition. 

4. RESULTS AND DISCUSSION 

4.1. Short timescale variability 

Except for two sources which are foreground Galactic 
stars, no source presents variability above the 99% confi- 
dence level. In Fig. 1 we show the results of the Bayesian 
block analysis for the only two sources that are variable 
on short timcscalcs. The dotted histograms represent the 
count rate binned in 400 second intervals. Note that the 
binning is for plotting purposes only. The green line is 
the power spectrum of the lightcurve. The thick solid 
line is the result of the Bayesian block analysis. The 
small square to the right of the lightcurve is a thumbnail 
image of the source from the observation. 

Source CXO J013327.7+304645 shows a clear X- 
ray outburst in ObsId 2023. This outburst, lasting about 
10,000 seconds completely accounts for the long term 
variability of the source. The lightcurve of the outburst is 
well fit by an exponential decay with a decay timescale of 
3900 seconds. The persistent luminosity in ObsId 2023 
is consistent within the errors with the luminosities in 
the other observations. Outburst and persistent counts 
arc not large enough to investigate changes in spectral 
shape. An X-ray color comparison of the burst and per- 
sistent emission shows a softening of the emission in HR2 
and a slight hardening in HRl that is significant only at 
the 2-sigma level. The optical counterpart to source CXO 
J013327.7+304645 is a star in the USNO catalog with a 
V magnitude of 16.8. We have analyzed archival HST 
WFPC2 data. The object was detected in bands cor- 
responding to U, B, V, and I. The colors are consistent 
with an early M type star, either on the main sequence 
or a giant. Using normal V magnitudes for M stars this 
puts the star at a distance of 90-320 pc for a main se- 
quence star, or 24-29 kpc for a giant. Considering the 
relatively high galactic latitude of M33, b = —31 deg, 
the lower distance value, and thus a main sequence star 
is more likely. Assuming a main sequence star the peak 
luminosity is between ^ 10^^ — 1.5 x 10"^° ergs s^^, and 
the total energy between ^ 10'^^ — 1.5 x lO '' ergs. 

Source CXO J013341.8+303848, a foreground star 
as well, has two outbursts in ObsId 786. The first 
outburst lasts about 2600 seconds, the second outburst 
18000 seconds later lasts about 800 seconds. The per- 
sistent level of X-ray emission increases slightly after 



each outburst. As shown in the long term lightcurve 
the source also has a long term trend to increasing lumi- 
nosity for the 3 observations. The decay of the outbursts 
cannot be fit uniquely; an exponential decay with a decay 
constant of ~1000 seconds or a linear decay are both pos- 
sible. During both outbursts the source becomes softer, 
but only at the 1 — 2a level. This source has no HST 
coverage. 

4.2. Long term variability 

Of a total of 261 sources, 198 have been detected in at 
least two observations, and 62 in all three observations. 
The luminosities for comparison are taken from spectral 
fits for the brighter sources (more than ~100 counts), or 
from the assumed spectrum, absorbed power law with a 
photon index of 2 and Galactic absorption. We find that 
49 of 198 sources show variability between observations 
above 3 sigma. An additional 29 sources show signs of 
variability between 2 and 3 sigma. Of the 49 variable 
sources only 16 are apparently persistent. 

The strength of the variability can in fact be used to 
confirm the nature of most of these sources as X-ray 
binaries. AGN in general show variability on weeks or 
months timescale only up to factors of 2-3 (Mushotzky 
et al. 1993; PaofiUo et al. 2004). Most of the sources in 
M33 show stronger variability, particularly since the flux 
ratios for most sources arc only lower limits. Moreover, 
the sources with flux ratios below 2 are bright sources 
that are unlikely to be AGN based only on brightness. 

Taking into account that ~80-90 of the 198 sources 
are likely background AGN, somewhat less than half the 
sources in M33 arc variable. Of these two thirds (34/49) 
may be candidate transients. This number is somewhat 
higher than the about 50% in the Milky Way and Mag- 
ellanic Clouds fLiu et al. 2000, 2001). But the detection 
limits (few 10^ erg s~^) do not allow to establish true 
transient behavior for undetected sources so the number 
of two thirds is only an upper limit. 

For the brighter sources {Lx > 5 x 10"^^ erg s~^) past 
X-ray satellites also provide data for long term variabil- 
ity. Fig. 2 shows data from Einstein HRI and IPC 
(Trinchieri et al. 1988), ROSAT HRI (Schulman & Breg- 
man 1995) and PSPC (Long et al. 1996), BeppoSAX 
(Parmar et al. 2001) as well as from this work for the 
bright sources in the fleld of view of Chandra . The 
quoted luminosities are converted to the Chandra band 
of 0.3-8.0 keV and, if they arc not individual measure- 
ments of spectra, converted to the spectral shape of a 
power law with photon index of 2 and the Galactic ab- 
sorption value towards M33, 6 x 10^° cm~^. The error 
bars in the plot only contain errors due to counting statis- 
tics. Other errors, e.g. due to the conversion of energy 
bands and different assumed spectral shapes, add gen- 
erally another 20-30% uncertainty in the luminosity. In 
addition different instrument responses, cross-calibration 
issues and other systematic effects add another source of 
systematic errors, which is hard to quantify. An excep- 
tion is the supernova remnant M33 X-14, for which con- 
verting the luminosity given by Long et al. (1996) would 
result in variability at over 6ct compared to other obser- 
vations. However, Long et al. (1996) assume a power 
law spectrum with a photon index of 2 to compute the 
source luminosity. Although this is a reasonable choice 
for X-ray binaries or AGN, the spectrum of M33 X-14 



4 



Grimm et al. 



01 3327.96+304644,8_2023,fits 
100 1000 10' 



01 3341 .88+303848.9_786.fits 
100 1000 





time (s) 



time (s) 



Fig. 1. — Left: Lightcurve and Bayesian block division of the outburst of CXO J013327. 7+304645, an X-ray active star. The panel 
shows the lightcurve binned in 400 second intervals (dashed histogram). Note that the lightcurve binning is for plotting purposes only. 
The Bayesian blocks are shown as the thick black histogram. The right part of the panel shows a thumbnail picture of the sonrce. Right: 
Lightcurve and Bayesian block division of the outburst of CXO J013341. 8+303848, another X-ray active star. 
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is quite soft. We extracted a spectrum of X-14 from 
the ROSAT observation rp600023a00, the longest of the 
PSPC observations with ~29 ks. The spectrum is well fit 
by an absorbed blackbody model with a column density 
of 4 X 10^^ cm~^ and a temperature of 0.09 keV. Al- 
though this model is rather unphysical for a SNR we are 
interested only in the flux, which is sufficiently accurate 
for the purpose of comparison. The ROSAT data also 
agree with the spectrum obtained from Chandra . Ex- 
trapolating this spectrum to the Chandra energy range 
gives a luminosity of 3.7 x 10^^ ergs s~^ which agrees very 
well with other observations of the source, in particular 
the Chandra observed value of 3.6 x 10"^® ergs s""'^. 

It is clear from Fig. 2 that four of the bright X-ray 
sources are variable. The sources that do not show ev- 
idence of variability at the 3cr level are M33 X-1, X-2. 
and the SNR X-14. 

The variability pattern is very similar to that observed 
in other X-ray source populations, regardless of galaxy 
type, e.g. The Antennae (Zezas et al. 2006), MlOl (Jenk- 
ins et al. 2005), NGC 4697 (Sivakoff et al. 2005). In 
particular that only very few, if any, sources show short- 
term variability is very common, although this is most 
likely due to limited photon statistics. On the other hand 
a large fraction (~10-40%) of sources exhibit long-term 
variability. Moreover, the fraction of sources with strong 
long term variability is quite confidently identified as X- 
ray binaries, and not background AGN. 

4.3. Spectra 

Fitting all the 43 sources detected with more than 
100 net counts with a power law and a thermal 
bremsstrahlung spectrum, with column density fixed 
to the Galactic value or as a free fit parameter, pro- 
vides about half of the sources with a good fit to the 
data. The other sources require either different mod- 
els, e.g. black body or other thermal plasma models, or 
multiple components. The results of the simple power 
law/bremsstrahlung fits support the validity of our as- 
sumptions of a fixed spectral model for conversion from 
counts to fluxes for fainter sources, assuming that the 
spectral properties of bright and faint sources are not 
systematically different (Grimm et al. 2005). 

In Fig. 3 we show histograms for the best fit power 
law and bremsstrahlung values for all sources. The up- 
per left panel shows the comparison between photon in- 
dices r of a power law with column density being a free 
fit parameter versus column density fixed to the Galac- 
tic value. The Galactic absorption towards M33 is only 
6 X 10^° cm~^. The upper right panel shows the same 
for bremsstrahlung temperature kT. The lower left panel 
shows the histogram of T for the case of Nh as a free 
fit parameter. The power law slopes are concentrated 
around 2, the canonical value for X-ray binaries and 
AGN. The peak at F = 5 comprises all sources with 
photon index larger than or equal to 5. Individually, all 
these sources are fit well with a blackbody or plasma 
model with temperatures of 0.1 0.13 kcV (blackbody) 
or 0.2 0.3 kcV (plasma model). The lower right panel 
finally shows the distribution of bremsstrahlung temper- 
atures. Note that the influence of Nh is stronger in case 
of a power law than for a bremsstrahlung spectrum as is 
evident in the generally smaller deviations from the one- 
to-one correlation for the bremsstrahlung temperatures 



compared to the photon indices in the upper panels. 

As expected, fits with a fixed low column density pro- 
duce smaller photon indices or larger bremsstrahlung 
temperatures, respectively. However, changes of tem- 
perature and photon index between fits with fixed and 
free column density are not significantly larger than the 
errors in the majority of cases. Comparing the differ- 
ence between photon indices divided by the square root 
of the errors shows that ^80% of the sources have values 
of less than 3. For thermal bremsstrahlung model the 
corresponding value is 73%. Thus the assumption of a 
general power law spectrum with F = 2, and Galactic 
absorption is quite good (Grimm et al. 2005). 

Using the results of the more detailed spectral fitting 
of the 43 bright sources, we can confirm the results of 
the relatively blind spectral fitting of all sources with a 
power law model. The left panel of Fig. 4 shows the 
photon indices of spectra that are well fit by cither a sin- 
gle power law (open circles) or a combination of a power 
law and other components (filled triangles) versus lumi- 
nosity. Errors are 90% errors on the slope. There is a 
clear trend that photon indices in multi-component fits 
are larger than for single power law fits, and there is an 
indication that brighter sources are more likely well fit 
by multi-component fits than by single power laws. The 
reason is most likely that the single power law sources 
have lower counts and multiple components are not dis- 
tinguishable. We also fit the single power law sources 
with a disk black body model, XSPEC model diskbb, 
which has the same number of degrees of freedom than 
the power law model. Two thirds of the spectra are well 
fit by a disk black body model, the best fit temperatures 
being in the range from 1 to 3 keV, as expected for X- 
ray binaries (Tanaka 2001). The inner disk temperatures 
versus luminosity are shown in the middle panel of Fig. 
4. The values are the same as for the sources that require 
disk black bodies for a good fit. This suggests that the 
single power law sources are likely to be X-ray binaries 
and only the low number of counts allows a good fit with 
a single power law. 

There are 10 spectra (5 sources) that are best fit 
with a single bremsstrahlung model with temperatures 
of '-^^0.1-0.3 keV. The bremsstrahhmg temperatures ver- 
s\is luminosity are shown in the right panel of Fig. 4 
(open squares). Since these values are unusually low for 
real bremsstrahlung sources, and the luminosities of the 
sources relatively low, 1 — 5 x 10^^ erg s"-*^, we fit these 
spectra with a thermal plasma model (XSPEC model 
a,pec) as well. With one exception we obtain good fits 
with the plasma model as well, the temperatures are in 
the range from 0.1 keV to 0.4 keV, and metallicities from 
0.05 to 0.3 solar with considerable uncertainties. These 
values are at or below the expected metallicity for M33 
(Blair & Kirshner 1985). The good fit quality is not sur- 
prising considering that the apec model has one degree 
of freedom more than the simple bremsstrahlung model. 
However, given the low temperatures and the consistent 
metallicities we consider the apec model to be the more 
physical model for these sources. The sources with good 
fits for the apec model are marked in Table Bl (Appendix 
B) and the apec model parameters arc given, in Fig. Bl 
(Appendix B) the contour plots for these fits are shown 
as well. Also shown in the right panel of Fig. 4 are 
six sources well fit by a black body model with very low 
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2. — Long term lightcurves of bright X-ray sources in M33 observed with Einstein, ROSAT, BcppoSAX, and Chandra . 



temperatures at or below 0.1 keV (filled triangles). For 
two sources the black body is the only component in 
the spectrum. These sources are candidates for super- 
soft sources and are discussed below in Sec. 4.4. Three 
of the other four sources have an additional power law 
component, that ranges from hard (F ~ 1.3) to very soft 
(F ~ 4.6). The black body temperatures are even lower 
than the temperatures inferred from ULX intermediate 
mass black hole candidates, see e.g. Miller et al. (2004). 
A truncated disk would be a possibility to explain the 
low temperature, similar to the scenario suggested by 
Kubota & Done (2004). However, given the luminosities 
the sources would be in the hard state (e.g. Maccarone 
(2003)) but, except for CXO J013444.6-I-305535, the pho- 
ton indices arc larger than 2. CXO J013444.6+305535 
could indeed be in the hard state with a photon index of 
1.3 but that number is not well constrained. This hard 
source could also be a candidate for a magnetic CV that 
are known to have hard spectra and relatively high lumi- 
nosities (Kuulkers et al. 2006) . The soft power law in the 
other sources might on the other hand be another ther- 
mal component. Alternatively, the soft emission might 
be generated in an outflow from the system, or heating 
of a surrounding medium. 

Fig. 5 shows a comparison of measured values for col- 
mnn density overlaid on a 1.49 GHz contour map of the 
central part of M33 from VLA (Condon 1987). The res- 
olution of the VLA map is about 1 arcminute, and the 



confusion limit is given as 0.1 mJy. There is no spatial 
correlation between the contour map and the magnitude 
of the measure colimin densities. Thus the X-ray absorp- 
tion value is in parts due to location of sources in front 
of/behind HI gas, and in parts due to intrinsic absorption 
around the X-ray source. 

To compare actual values for the column density due 
to HI, we compute the brightness temperature of the HI 
gas according to 



_2fct^ 



(2) 



with Si, the radio flux in units of Jy/steradian at fre- 
quency V (1.49 GHz), fl the opening angle of the beam 
(1'), and Tb the brightness temperature; k is the Boltz- 
mann constant, c the speed of light. Assuming optically 
thin emission we compute the HI column density accord- 
ing to 

Nh = 1.83 X 10^^Tb[Kcm-% (3) 

Within the size of a beam the radio flux has an RMS of 

^5 mJy around the zero point. Therefore we take 5 mJy 
as the upper limit on the sensitivity in a beam. This 
flux corresponds to ^320 K or Nh of 5.7 x lO'^" cm^^, 
which is very close to the Galactic absorption value. Fig. 
6 shows the comparison between the X-ray absorption 
values from spectral fitting and the "expected" absorp- 
tion from HI. Filled circles are sources in M33, crosses 
are a selection of Galactic X-ray binaries from Vrtilek 
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Fig. 3. Overall spectral properties of the samijle of 43 bright X-ray sources in M33. The panels show the correlations of photon index 
r with Nf{ free versus fixed (upper left), breinsstrahluiig temperature kT with free versus fixed (upper right). Note that the influence 
of Nh is stronger in case of a power law than for a bremsstrahlung spectrum as evident in the smaller deviations from the one-to-one 
correlation for the bremsstrahlung temperature. The lower panels show the distributions of F (lower left) and kT (lower right) for the case 
of Nff as a free fit parameter. The power law slopes are mainly in the range from 1.4 to 2.5, the canonical values for X-ray binaries and 
AGN. 



et al. (1991). Upper limits are shown as arrows. The 
vertical and horizontal lines denote the Galactic absorp- 
tion value. The figure can be separated in roughly three 
regions: 1. Sources at or below the horizontal line. i.e. 
sources with apparently less absorption than the Galactic 
value, which probably have an undetected thermal com- 
ponent in the spectrum. 2. Sources below the diagonal 
line which are probably located in front of at least parts 
of the HI emission. These sources arc most likely on the 
Milky Way facing side of M33. 3. Sources above the di- 



agonal line are intrinsically absorbed. No conclusion as 
to their location can be drawn. 
Values for Galactic X-ray binaries are taken from Vr- 

tilek et al. (1991) and Liu ct al. (2001) for Ay values 
missing in the former. We convert the optical extinction 
to neutral hydrogen column density according to Predehl 
& Schmitt (1995) A^^ = 1.77 x Ay in units of lO^^ cm'^, 
and use Ay = 3 x E{B — V) in cases where only the red- 
dening in available. The comparison with Galactic X-ray 
binaries shows that M33 X-ray sources also often show 
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Fig. 4. — The left panel shows the best fit photon indices for sources that are well fit by a single power law (open circles) or a multi- 
component spectrum containing a power law (open circles with filled triangles) versus luminosity. Two thirds of the single power law 
sources are also well fit by a disk black body model. The middle panel shows the disk black body temperatures versus luminosity for 
sources that are well fit by a single power law, but also by a simple disk black body model (filled circles). As expected for X-ray binaries 
the disk temperatures are in the range from 1-3 keV. The temperatures are also in the same range as the temperatures for sources that 
require a disk black body component in the spectrum (open triangles). The right panel shows temperatures of sources that are either well 
fit by a bremsstrahlung model (open squares) or a black body model (filled triangles) versus luminosity. With the exception of two spectra, 
the bremsstrahlung sources are also well fit by a thermal plasma model, XSPEC model apec, as discussed in the text. Errors are 90% 
confidence level. 



excess absorption at levels comparable with Galactic X- 
ray binaries. The larger values for A''// (HI) in the Milky 
Way can be explained simply by our location in the Milky 
Way, compared to the nearly face-on orientation of M33. 

Only six of the 25 sources observed in more than one 
observation with more than 100 counts show significant 
variability in the hardness ratio diagram, see Appendix 
B. Three of the sources are the well known sources, M33 
X-4, M33 X-5, and M33 X-13. The other three are CXO 
J013329.0-h304216, CXO J013329.3-H304508, and CXO 
J013410. 5-1-303946. There is a clear separation in that 
most sources have a very small significance for moving 
in the hardness ratio diagram, whereas four of the ones 
that show changes in hardness ratio are above 99.99% 
confidence level. 

Based on spectral shape (photon indices, tempera- 
tures) and absorption properties we conclude that M33 
X-ray sources as a population are very similar to Galactic 
X-ray binaries. 

4.4. Individual sources 

Here we discuss the properties of individual bright M33 
sources. These include previously known sources (M33 
X-...) and sources newly discovered in our Chandra ob- 
servations. Wide band spectral studies exist in the litera- 
ture for some of these sources. In particular, BeppoSAX 
has observed the brightest sources in M33 and measured 
spectra of the sources M33 X-4 to M33 X-10 in the en- 
ergy range from 2.0-10.0 keV (Parmar et al. 2001). Of 
these sources X-6 and X-10 are not in the field of view of 
the Chandra observations, and X-8 is the nucleus which 
is not discussed here. 

M33 X-4 (CXO J013315. 1+305317): The X-ray 
spectrum of this soft source measured by BeppoSAX is 
well fit by either a power law model with photon index 
of ~3 or a bremsstrahlung model with kT « 2 keV. This 
source has been observed by Chandra twice, in Obslds 
786 and 2023. Between the observations, the flux in- 
creased by a factor of two from 5 x 10'^^ erg s~^ to lO'^* 
erg s~^. The variability rules out the interpretation of 
the source as a SNR (Okada et al. 2001). Using the 



overlapping energy range with the BeppoSAX observa- 
tions, 2.0-8.0 keV, we obtain the same values for absorp- 
tion, photon index, and bremsstrahlung temperature well 
within the errors and x^-values below 1 for observation 
786. For observation 2023, although the best fit values 
are consistent with the BeppoSAX and Obsid 786 values, 
the x^-values are well above 2. Moreover, for the whole 
Chandra range, 0.3-8.0 keV, none of these models is a 
good fit. A statistically satisfactory fit for both Obslds 
can be achieved with a model consisting of a black body 
and a power law. The black body temperatures ('^0.65 
keV) and photon indices (3-4.7) for the two observations 
are consistent within the errors. However, the absorption 
to the source changes. Taking into account the different 
values for absorption there is a clear distinction at 99% 
confidence level in the contour plots, for T versus Nh and 
kT versus Nh ■ The difference in F versus kT is significant 
only at 90% confidence level. The black body tempera- 
ture is relatively low, and does not change significantly, 
but similar values for disk black body temperatures have 
been observed in Galactic BHs. These temperatures oc- 
cur at luminosities of ^^1-5% of the Eddington luminosity 
(Gierliiiski & Done 2004). With a luminosity increasing 
from ~ 5 X 10^^ erg s~^ to 10^^ erg s~^, this interpreta- 
tion would make the compact object in X-4 a rather mas- 
sive but still stellar mass BH given the luminosity. The 
photon index of the power law component is much softer 
than expected for a low-hard state black hole (Tanaka 
2001) but is in agreement with high-soft state Galactic 
black holes (McChntock & Remillard 2006). The change 
in the photon index between the observations (from ^3 
to 4.7) does not indicate a state change from low-hard 
to high-soft. Even a photon index of 3 is too large for 
the low-hard state (F ~1.7). The photon index is also in 
the range observed for the Galactic XRB Sco X-1 which 
has a neutron star primary (Bradshaw et al. 2003), but 
Sco X-1 does not exhibit a soft component like M33 X-4 
(Kahn et al. 1984). Thus based on the X-ray data the 
source is likely to be a BH binary. 
The hardness ratio analysis shows that the hardness of 
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Fig. 5. — Measured values for column density overlaid on a 1.49 GHz contour map of the central part of M33 from VLA (Condon 1987). 
The resolution of the VLA map is about 1 arcminute, and the RMS is about 0.1 mjy. Prominently visible is the HII region NGC 604 at 
the upper left. Also visible towards south of the center is the southern spiral arm. There is no correlation between the contour map and 
the magnitude of the measure column densities. 



the source varies significantly between the observations. 
The source becomes harder in HRl and softer in HR2, 
while the luminosity increases by a factor of ~2. This 
change is due to the increase in the medium band (1.0- 
2.1 keV) with respect to the soft and hard band, which 
is reflected in the spectral fits by a higher absorption 
column and a lower power law photon index. 

M33 X-5 (CXO J013324.4+304401): This source 
has been observed in all three Chandra observations. 
For M33 X-5 (CXO J013324.4+304401) BeppoSAX as 



well finds good fits with a soft power law (F « 3) or 
bremsstrahlung model {kT « 2.8 keV). In all Chan- 
dra observations the overlapping energy range is fit sat- 
isfactorily with these values. However, for the whole 
Chandra energy range neither a single power law nor a 
bremsstrahlung spectrum give good fits. A good fit is ob- 
tained with a disk black body (~1.3 keV) and power law 
(F = 4 — 9) combination. The contour plots in Fig. ?? 
(Appendix B) show no significant variation of the spec- 
tral parameters given the large uncertainties. The er- 
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Fig. 6. — Measured values for column density from X-ray spectral fitting versus absorption values inferred from HI radio flux. Upper 
limits are shown as arrows. Vertical and horizontal lines denote the Galactic absorption value. Filled circles are sources in M33, crosses are 
Galactic X-ray binaries, see text. The plot can be separated in three regions: 1. At/below the horizontal line, sources with apparently less 
absorption than the Galactic value. 2. Below the diagonal line, sources are probably located in front of at least parts of the HI emission. 
3. Above the diagonal line, sources are also intrinsically absorbed. 

rors on the photon index are large enough to not make 
the changes between observations significant. However, 
a disk black body alone does not yield a good fit. The 
disk black body temperature is in good agreement with 
observations of Galactic BH binaries at moderately sub- 
Eddington luminosities (Gierlihski & Done 2004). The 
luminosity of the sources changes between 3 x lO'^^ erg 
s~^ and 5 x 10^^ erg s~^; thus it is significantly below 
the Eddington limit for a black hole but in the range of 
transition luminosities from the low-hard to the high-soft 
state of Galactic black holes (Maccarone 2003). Another 



interpretation is that the source is a neutron star and 
the steep power law corresponds to the harder thermal 
component from the neutron star surface (Tanaka 2001). 
Due to the complexity of the combined model only the 
disk temperature is well determined. 

However, the hardness ratio analysis shows significant 
variation between Obslds 1730 and 786 on the one hand 
and Obsid 2023 on the other. X-5 becomes significantly 
softer in HRl but does not change significantly in HR2. 
There seems to be no correlation with luminosity as Ob- 
slds 786 and 2023 have the same luminosity but signifi- 
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cantly different HRl values. 

M33 X-7 (CXO J013334.1+303210): This source 
has been observed only in Chandra observation 1730. 
For M33 X-7 (CXO 013334.1-1-303210) BeppoSAX as 
well finds good fits with a power law (F « 2 — 3) 
or bremsstrahlung model {kT « 3.7 keV). In the 2.0- 
8.0 kcV energy range both models provide satisfactory 
(power law) to good (bremsstrahlung) fits to the Chan- 
dra data. The photon indices for fixed and variable 
absorption are larger than in the BeppoSAX obser- 
vations but still within the rather large errors. The 
bremsstrahlung temperatures for fixed and variable ab- 
sorption arc smaller, and in the case of fixed absorption 
inconsistent with the lower limit of the BeppoSAX obser- 
vation. Although the errors are quite large this indicates 
that the source was softer in the Chandra observations. 
Moreover, M33 X-7 was about five times brighter in the 
Chandra observation (~ 7 x 10'^'' erg s~^) compared to 
the time BeppoSAX observed the source (~ 1.6 x 10^^ 
erg s~^) (note that both luminosities are based on spec- 
tral fits). The only satisfactory fit to the whole energy 
range of Chandra is an absorbed disk blackbody with 
an inner disk temperature of 1 ± 0.02 kcV. This value is 
consistent with inner disk temperatures of Galactic black 
holes (McChntock & Remillard 2006). The column den- 
sity is not well constrained, but below 6.5 x 10^° cm~^, 
putting the upper limit of the absorption at the Galac- 
tic value. Pietsch et al. (2004) obtain the same values 
within the errors for an XMM observation, and also new 
Chandra observations give the same spectral parameters 
(Pietsch et al. 2006). This might indicate the presence 
of an additional soft component that due to insufficient 
counts presents itself as a low absorption value. An ad- 
ditional soft component, however, does not improve the 
fit, and its parameters are not well determined. Fig. 7 
shows the spectrum and the contour plot of absorption 
column density versus inner disk temperature. Given 
the softer spectrum in the Chandra observation, a possi- 
ble additional soft component, and the five times higher 
luminosity, it is possible that X-7 underwent a state tran- 
sition. On the other hand the photon index of the power 
law in the BeppoSAX spectrum is already softer than 
expected for a low-hard state source (Tanaka 2001). 

CXO J013343.4+304630 and CXO 

J013409.9-I-303219: These sources are the best 
candidates for super-soft sources (SSS) in M33. CXO 
J013343.4+304630 is well fit by a black body spectrum 
with a temperature of about 75 eV and moderate 
absorption (Njf 8 x 10^^ cm^^). The source shows 
no significant spectral variability, but is variable on 
timescales of months. CXO J013409.9-F303219 has been 
observed only once. Its spectrum is also well fit by 
a black body with a temperature of 37 eV and also 
moderate absorption (Nh ^ 7 x 10^^ cm^^). The tem- 
peratures are in the range for SSSs, and the sources are 
not significantly detected above 1 keV, see e.g. Kahabka 
& van den Heuvel (1997). Given their luminosities 
of ~ 3 - 10 X 10^5 erg s'^ (CXO J013343.4-h304630) 
and - 3 X 10^6 erg s^i (CXO J013409.9+30321) these 
sources are most likely nuclear burning white dwarfs. 

M33 X-9 (CXO J013358.8+305004): This source 
was observed in Chandra observations 786 and 2023. It 
is well fit by a power law with photon index of ~1.2 in 
the BeppoSAX observation. For a bremsstrahlung model 



only lower limits to the temperature are provided (> 2.6 
keV). It is important to note that in the BeppoSAX ob- 
servations X-9 is unresolved, but ROSAT observations 
have shown that X-9 actually consists of 3 sources, at 
least two of which seem to be variable and of compara- 
ble flux (Long et al. 1996). Since BeppoSAX does not 
resolve the sources, it is impossible to say if any of them 
dominate the spectrum and if yes which one. The Chan- 
dra source is associated with the source X-9a based on 
the ROSAT position. In the whole Chandra energy range 
a power law is a good fit. The absorption column is not 
well constrained in this fit, but the values for the pho- 
ton index are ~ 1.7, with a reduced of 0.87. Fig. 8 
shows the spectrum and the contour plot of absorption 
column density versus photon index. Assuming for the 
sake of the argument that a significant fraction of the 
BeppoSAX fiux originated from X-9a, the source was 
thus significantly harder during the BeppoSAX obser- 
vation. Even if there is no spectral change in X-9 the 
source was then about a factor of ten brighter during the 
ROSAT observations where it reached ~ 1.4 x 10^'^ erg 
s""'^. Between the two Chandra observations the source 
luminosity decreases by a factor of ^2 from 1.7 x 10^^ 
erg s~^ to 7 X 10^^ erg s~^ without a significant change 
in spectral parameters. Based on the variability, either 
spectral and/or timing, and the luminosity the source is 
thus highly likely to be an X-ray binary. 

M33 X-13 (CXO J013354.8+303309) and CXO 
J013329. 0-1-304216: These two sources have very soft 
bremsstrahlung temperatures (~ 0.2 keV) but are also 
variable on timescales of months, thus excluding SNRs 
or HII regions as counterparts. 

M33 X-13 (CXO J013354.8+303309) is weU fit with a 
bremsstrahlung model with temperature 0.22 keV and 
absorption of ~ 3 x 10^^ cm~^ in both observations. The 
source is coincident with a SNR (Gordon et al. 1999), 
but the fiux varies by a factor of two between two obser- 
vations. It also shows significant change in the hardness 
ratios. The source becomes softer in HRl and harder in 
HR2. The flattening of the spectrum is accompanied by 
an increase in luminosity 2 x 10^^ erg s~^ to 4 x 10^^ erg 
s-i. 

CXO J013329.0-f 304216 shows a significant softening 
in HRl and hardening in HR2 from Obsid 1730 to Ob- 
sld 2023, the same behavior at similar hardness ratios 
like M33 X-13 (CXO J013354.8-I-303309). The spec- 
trum of the source is well fit by a bremsstrahlung spec- 
trum with a temperature decreasing from 0.24±0.01 kcV 
(ObsId 1730) to 0.21 ± 0.006 kcV (ObsId 786) and to 
0.19 ± 0.006 keV (ObsId 2023). However, the temper- 
ature is somewhat degenerate with the column density 
which increases from (0.29 ± 0.1) x 10~^^ cm~^ (ObsId 
1730) to (0.36 ± 0.1) x 10-22 cm^^ (Qbsld 786) and to 
(0.42 ± 0.17) X 10-22 cm-2 (ObsId 2023). Although the 
column density changes systematically, the large errors 
make it impossible to choose which, if any, value is cor- 
rect. Thus although the errors on the temperatures are 
small the significance of the temperature change from the 
spectral fit is not large. Between the Chandra observa- 
tions the luminosity increases from 2.7 x 10'^^ erg to 
3.6 X lO^'^ erg s'^ 

Temperature and variability of these sources indicate 
that they could be quasi-soft sources, although a sim- 
ple black body does not provide a good fit. However, 
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Fig. 8. — Energy spectrum and contour plot for M33 X-9. The model spectrum is an absorbed power law with a photon index of 
1.88 ± 0.04. The reduced is 0.87. 




more complicated models (e.g. black body plus power 
law) do not improve the fits due to the small number of 
observed counts. For quasi-soft sources the luminosities 
would be relatively low (a few 10^^ erg s~^) but still in 
the range observed in other galaxies (Di Stcfano & Kong 
2004). The sources are too bright 2 - 4 x lO^^ erg s"^) 
for normal CVs xo^°~^^ erg s~^), and magnetic CVs 
that can reach such luminosities have very hard spec- 
tra with effective temperatures of several keV (Kuulkers 
et al. 2006). On the other hand the sources have too high 
temperatures for super-soft sources (Kahabka & van den 
Heuvel 1997). 

CXO J013329. 3+304508 and CXO 

J013410.5+303946: The two sources show evi- 
dence for spectral variability in their hardness ratio 
diagrams. Source CXO J013329.3-I-304508 shows a 
significant softening in HRl but no change in HR2 
from Obsid 1730 to Obsid 2023. The best fit model for 
the spectrum is a bremsstrahlung model, however the 
temperature is not well determined in the fit. Moreover, 



the source is variable between Obslds 786 and 2023, 
between which the luminosity decreases by ~30% from 
3 X 10^6 erg s'^ to 2 x lO^^ erg 8"^ Source CXO 
J013410.5-f 303946 shows a softening in HRl from Obsid 
1730 to Obsid 786 but no significant change in IIR2. 
From Obsid 786 to Obsid 2023 the HRl hardens again 
so that the hardness ratio for Obsid 2023 is consistent 
with the hardness ratio for Obsid 1730. The spectrum 
of Obslds 1730 and 2023 is well fit by a disk black body 
with an inner disk temperature 1.2 keV. The spectrum 
of Obsid 786 on the other hand is well fit by a power law 
with a photon index of 1.7. Both photon index and disk 
temperatures are consistent with Galactic X-ray binary 
spectra (Tanaka 2001). It is interesting to note that the 
luminosity increases between the first two observations 
and stays high in the third, although the spectral shape 
in the first and third observation are basically identical. 
The source increased its luminosity between Obslds 1730 
and 786 from ~ 1.3 x 10^^ erg s'^ to ~ 2.1 x 10^'^ erg s'^ 
and had the same high luminosity during Obsid 2023. 
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Although the power law spectrum is relatively hard, 
the softening of the hardness ratio and the increasing 
luminosity do not indicate a state change in the system. 

The spectral and timing behavior of the bright X- 
ray sources in M33 is varied but can be well described 
within the framework of Galactic X-ray sources and X- 
ray sources in other galaxies observed with Chandra . 

5. SUMMARY AND CONCLUSION 

In the three Chandra observations of M33, performed 
between the years 2000 and 2001, 261 sources have 
been detected in the range of X-ray luminosities from 
Lx ~ 1034-38 erg s^^ (0.3-8.0 keV). Of this total, 198 
sources have been detected in at least two observations, 
and 62 sources have been observed in all three observa- 
tions. We find that 49 sources show variability between 
observations above 3 sigma. For a total of 43 sources, 
25 of which in more than one observation, the number 
of counts is sufficient for a more detailed spectral fitting. 
Given the angular extent of our survey we expect ~3- 
4 AGN interlopers for sources with luminosities above 
IQ^e gj-g g-i^ g^^^ about one for sources with luminosities 
above 10^^ erg s~^ (see Grimm et al. (2005)). All sources 
in M33, except for the nucleus, have luminosities below 
~ 10^8 erg s~^ 

Time variability: Roughly a quarter (49/198) of 
all Chandra soiirces exhibits long timescalc variability 
(^wecks years), therefore we can exclude SNRs and HII 
regions in M33 as their counterparts. Except for two 
Galactic star interlopers none of the M33 sources show 
any short timescale X-ray variability (~seconds-hours) . 
Detection of such variability would separate stellar mass 
objects from AGNs, since in AGNs the light travel time 
is too large to produce variability on time scales shorter 
than hours. However, the amplitude of the long-term 
variability (larger than a factor of 2-3) makes it unlikely 
for the variable sources to be AGN (generally less than 
a factor of 2-3). 

A significant number of sources showing long-term 
variability but little to no short-term variability is very 
similar to other extragalactic X-ray source populations. 
The comparison of the amplitude of long term variability 
with AGN and known X-ray binaries shows that a large 
majority of these sources arc X-ray binaries. 

X-ray spectra: The large majority of source spec- 
tra is well fit by models containing either a power law 
or a thermal plasma (bremsstrahlung or apec), but in 
the most luminous sources more complex models are re- 
quired. Of the 43 sources with enough counts for detailed 
spectral modeling, 

• ten are well fit with spectral models consistent 
with accreting X-ray binaries, namely a disk black 
body or a black body plus power law. Three of 
the sources well described with a black body plus 
power law combination have black body temper- 
atures at or below 0.1 kcV. These are even lower 
than temperatures inferred from ULX interinediate 
mass black hole candidates, see e.g. Miller et al. 
(2004). The luminosities of these sources also do 
not indicate the presence of an intermediate mass 
black hole. 

• In addition, 16 sources are well fit by a single power 
law. With one exception all the photon indices are 



in the range ^1-2, consistent with Galactic X-ray 
binaries (Tanaka 2001). The photon indices are 
not inconsistent with AGN but we expect only ~3 
background AGN among the bright sources. Two 
thirds of the sources well fit by a single power law 
are also well fit by a disk black body model with 
inner disk temperatures in the range 1.0 3.0 keV, 
which would be consistent with Galactic BH and 
NS binaries (Tanaka 2001). 

• There are also 4 sources with bremsstrahlung spec- 
tra and temperatures between ^2-7 keV. 

Moreover, 12 of 30 sources also show long term variabil- 
ity, confirming their nature as accreting objects. There- 
fore more than half (30) of the sources are consistent with 
being X-ray binaries based on the X-ray spectrum. 

The majority (12) of the remaining 13 sources are 
well fit by thermal plasma models (bremsstrahlung or 
apec models as discussed above), 7 of which with rela- 
tively low temperatures around 0.2 keV. Five of these 12 
sources are variable. Three sources with unusually low 
bremsstrahlung temperatures from 0.1-0.3 keV are also 
well fit by an apec model as well; one is variable. The 
apec temperatures of these sources are around 0.3 keV 
and abundances are consistent with the overall metallic- 
ity of M33 (Blair & Kirshner 1985) (see Fig. 4). A black 
body model docs not provide a good fit to the data. Five 
non-variable of the 12 sources are coincident with SNRs 
(Gordon et al. 1999). Variable sources in this category 
are candidates for quasi-soft sources. 

Two sources are well fit by a single black body with 
temperatures of ^75 eV and ^37 eV respectively. One 
of the sources is a long term variable with no significant 
spectral variability. The other source has been observed 
only once. Both sources are not significantly detected 
above 1 keV. These sources are the best candidates for 
super- soft sources in M33. 

Given the spectra, and variability, the majority of these 
sources are likely X-ray binaries, with a sizable fraction 
of SNRs/HII regions confirming a young X-ray source 
population. 

Spectral variability: Six of the bright sources show 
evidence for spectral variability. The type of variability 
is varied, as has been seen in other galaxies as well (Zezas 
et al. 2006). The spectral variability is sometimes corre- 
lated with luminosity (M33 X-4, X-13), and sometimes 
not (M33 X-5). In the source X-4 and X-5 the spectral 
change seems to originate in the column density-photon 
index degeneracy. In both cases, the thermal compo- 
nent (black body or disk black body) remains constant, 
whereas the photon index and column density change 
their values. With the current data it is not possible 
to decide, if there is a change in column density within 
the system or a genuine change in the hard component, 
assuming that the model is a correct description. M33 
X-13 and CXO J013329.0-F304216 both show a continu- 
ous softening of the spectrum that is correlated in both 
cases with an increase in luminosity, and they could be 
identified with quasi-soft sources. The other two sources 
show significant change only in the soft band with only 
moderate change in luminosity. 

Given the brightness, spectra, and variability (tempo- 
ral and spectral) these sources are highly likely X-ray 
binaries. 
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A comparison of Nu values measured from the X-ray 
spectra and HI radio observations of M33 shows that the 
X-ray absorption is larger for most sources than inferred 
from HI column density. A few sources have lower X-ray 
absorption than expected and may be located in front 
of the HI gas. However, the majority of sources is in- 
trinsically absorbed. The absorption in X-rays can be 
compared with values observed in Galactic X-ray bina- 
ries, and is generally moderate at a few 10^^ cm~^. 

This comparison allows the localization of a few sources 
in M33, and indicates for most sources intrinsic absorp- 
tion at the same level as observed in Galactic X-ray bi- 
naries. 

The luminosity function of M33 sources and optical 
counterparts from Grimm et al. (2005) have shown that 
the X-ray source population in M33 is dominated by 
young objects. The long term variability, spectra, and 
spectral variability presented in this paper provide evi- 
dence that a large majority of sources are X-ray binaries. 
The properties of the X-ray binaries (variability, spec- 
tral parameters including absorption) are very similar 
to HMXBs in the Milky Way or the Magellanic Clouds. 
Thus X-ray binary populations do not show strong vari- 
ations over the parameters, e.g. metallicity, covered by 
the Milky Way, M33, and the Magellanic Clouds. Thus 
the M33 X-ray source population confirms our knowledge 
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about HMXB populations in the Milky Way and other 
galaxies. Therefore M33 can also be used as a comparison 
template for similar galaxies at larger distances, which is 
important as galaxies of this type are numerous. 

Comparisons with previous X-ray missions show that 
repeated observations of a galaxy are important to study 
the X-ray source population in detail. This analysis also 
shows that especially with repeat observations it is possi- 
ble to identify the major parts of an X-ray population on 
X-ray data alone. The ongoing deep Chandra survey of 
M33 will add significantly more data to such an analysis, 
and provide a more detailed picture of the M33 X-ray 
source population. 
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A. TIME VARIABILITY 

Table Al presents the analysis results of variability between observations. We use Poissonian statistics for the 
background-subtracted source counts. The table shows only sources with variability above 3a between any two ob- 
servation. The first column gives the name of the source, the next three columns give the fluxes, with errors in 
parentheses, for each observation if the source was detected, or 3 a upper limits on the flux. The last three columns 
give the significance of variability between any two observations. Comparisons between upper limits are set to zero. 

TABLE Al 

List of sources variable at more than 3<t between any observation. Comparisons between upper limits are set to zero. 
Source Flux [10~^ cts/s/cm"^] Significance [cr] 









1730 


786 


2023 


1730 
786 


786 
2023 


1730 
2023 


cxo 


J013308, 


.3+304802 


— 


66 (8.4) 


<4.5 


— 


7.3 


— 


cxo 


J013315, 


.1+305317 


— 


3000 (46) 


5800 (62) 


— 


36.3 


— 


cxo 


J013315, 


.5+304448 


180 (18) 


170 (19) 


<4.3 


0.4 


8.7 


9.8 


cxo 


J013321. 


.7+303858 


24 (5.2) 


— 


<6.2 


— 


— 


3.4 


cxo 


J013321, 


.9+303921 


22 (5.0) 


— 


<4.3 


— 


— 


3.5 


cxo 


J013323, 


,9+304821 


— 


46 (7.4) 


<9.2 


— 


5.0 


— 


cxo 


J013324, 


,4+304401 


1700 (39) 


2800 (45) 


2500 (31) 


18.5 


5.5 


16.1 


cxo 


J013327, 


,7+304645 


14 (4.3) 


20 (5.4) 


52 (6.6) 


0.9 


3.8 


4.8 


cxo 


J013328, 


,6+304321 


20 (4.7) 


19 (4.9) 


<3.5 


0.1 


3.2 


3.5 


cxo 


J013329, 


,0+304216 


150 (12) 


160 (11) 


200 (10) 


0.6 


2.7 


3.2 


cxo 


J013329, 


,2+304537 


34 (6.0) 


39 (6.5) 


<6.2 


0.6 


5.0 


4.6 


cxo 


J013329, 


,2+304508 


140 (11) 


170 (12) 


120 (8.2) 


1.8 


3.4 


1.5 


cxo 


J013333, 


,0+304920 


— 


28 


<6.8 


— 


3.5 


— 


cxo 


J013334, 


,1+303714 


43 (6.4) 


— 


<4.3 


— 


— 


6.0 


cxo 


J013335, 


,5+303728 


12 (4.0) 


19 (3.7) 


<4.3 


1.3 


4.0 


1.9 


cxo 


J013336, 


,3+303742 


90 (9.1) 


48 (5.5) 


<7.3 


3.9 


7.4 


9.1 


cxo 


J013337. 


,4+304718 


71 (8.4) 


100 (9.7) 


29 (5.2) 


2.3 


6.5 


4.3 


cxo 


J013339, 


,2+304049 


50 (6.7) 


47 (8.4) 


<6.2 


0.3 


4.9 


6.5 


cxo 


J013340, 


,0+304323 


32 (5.6) 


23 (4.7) 


<9.5 


1.2 


2.9 


4.0 


cxo 


J013341, 


,5+304136 


26 (5.0) 


37 (5.6) 


<10.7 


1.5 


4.7 


3.1 


cxo 


J013341, 


,8+303848 


70 


150 


200 


6.8 


3.8 


10.4 


cxo 


J013342, 


,5+304253 


15 (4.2) 


56 (6.8) 


<11.8 


5.1 


6.5 


0.8 


cxo 


J013343, 


,4+304630 


20 (8.8) 


61 (7.7) 


79 (7.5) 


3.5 


1.7 


5.1 


cxo 


J013344, 


,2+304026 


12 


30.8 


<5 


2.3 


3.9 


1.5 


cxo 


J013346. 


,2+303807 


<3.5 


<4.3 


90 (9.0) 


0.0 


9.5 


9.6 


cxo 


J013350. 


,5+303821 


27 (5.1) 


57 (5.7) 


<11.8 


3.9 


7.9 


3.0 


cxo 


J013353, 


,6+303605 


<10.4 


9.9 (2.9) 


— 


0.2 


— 


— 


cxo 


J013354, 


,8+303309 


120 (10) 


2.30 (12) 


— 


7.0 


— 


— 


cxo 


J013356, 


,8+303729 


37 (8.4) 


83 (6.7) 


— 


4.3 


— 


— 


cxo 


J013358, 


,8+305004 


— 


93 (9.8) 


36 (4.5) 


— 


5.3 


— 


cxo 


J013409, 


,8+305044 


— 


3.8 


33 (4.5) 


— 


5.3 


— 


cxo 


J013410, 


,3+305346 




<5 


23 (3.8) 




4.7 




cxo 


J013410, 


,5+303946 


72 (7.9) 


120 (8.0) 


120 (7.7) 


4.3 


0.0 


4.4 


cxo 


J013416, 


,7+305101 




<7.7 


27 (3.7) 




5.2 




cxo 


J013419, 


,2+304942 




<5 


17 (3.0) 




4.0 




cxo 


J013424, 


,6+304428 


<5.1 


26 (6.0) 


<6.8 


3.5 


3.2 


0.0 


cxo 


J013425, 


,3+304157 


<15.2 


33 (6.2) 


11 (2.9) 


2.9 


3.2 


1.4 


cxo 


J013426. 


,7+304811 




<3.5 


14 (2.8) 




3.8 




cxo 


J013429, 


,1+304212 


13 (4.5) 


<14.7 


7.5 (2.2) 


0.4 


3.3 


1.1 


cxo 


J013429, 


,7+305026 




<7.7 


34 (3.9) 




6.7 




cxo 


J013432, 


,1+305158 




<5 


93 (13) 




6.8 




cxo 


J013432, 


,5+305035 




<5 


15 (3.0) 




3.3 




cxo 


J013432, 


,7+303436 


51 (8.3) 


<3.6 




5.7 






cxo 


J013433, 


,7+304701 




430 (27) 


280 (14) 




4.9 




cxo 


J013435, 


,0+304439 




<7.7 


17 (3.0) 




3.1 




cxo 


J013436, 


,4+304713 




50.3 


23 (3.3) 




3.1 




cxo 


J013446, 


,7+304449 




<6.4 


21 (3.3) 




4.4 




cxo 


J013449, 


,0+304446 




<6.4 


34 (4.0) 




6.9 




cxo 


J013451, 


,1+304356 




<3.5 


30 (3.8) 




7.0 




cxo 


J013451, 


,9+304615 




120 (14) 


48 (4.8) 




4.9 
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B. SPECTRA 

The spectra are combined in one plot. Black corresponds to the first observation, red to a second observation, and 
green to a third obcrvation. For easier comparison contour plots arc also overlaid for all observations. Solid contours 
correspond to the first observation, dashed contours to the second observation, and dotted to the third observation. 
Note that not all sources are observed in all observations. 

Some sources are well-fit by a bremsstrahlung model with low temperatures. We also fit these spectra with an APEC 
model which also provides a good fit. In case the APEC model provides a good fit as well we give the parameters and 
contour plots in the table and figures. 



TABLE Bl 
Spectral parameters. 



Source 


Obsid 


XSPEC model 


Nh 
[10^2 cm-2] 


parameters 
kT, TO, Tin, E, a: [keV], 
Z: [solar abund.] 










kT 




CXO J013253.4+303817 


1730 


bremss 


"^^—0.03 


7 1 +3.55 
'■-^-2.0 




(M33 X-1) 


















i 




CXO J013253.9+303312 


1730 


pow 


16+°°^ 


1 Q+0.26 
-l-='-0.23 




(M3ci X-2) 














kT 




CXO J013308.3+304802 


786 


bremss 


17+0-12 


n 0+4.3 












kT 


r 


CXO J013315. 1+305317 


786 


bbody+pow 


0-34Io:i3 


0.62ini 


o n+0.3 
ci.U_Q 7 


(M33 X-4) 


2023 


bbody+pow 




0.69+0.04 


. „+0.75 


CXO J013315.5+304448 


1730 


pow 




r 

9 7+0.75 
^■'-0.65 






786 


pow 


43+1.02 


5 8+2-'' 












Tin 


r 




1730 


n 1 oH-ir^— 1— T>mnr 
UlOiiUUT^pUW 


< 0.56 


1 O+0.15 
^••^-0.11 


9.5" 


(M33 X-5) 


786 


diskbb+pow 




1 Q+0.23 
^■•^-0.1 


r. 7+1.3 
^•'-1.9 




2023 


diskbb+pow 




1 4+0- 13 

-L.I Q 2 


3.5+Q 1 










kT 


r 




2023 


DDOdy-rpow 


\ u.oo 


1.2" 


„ ^+1.7 










kT 




CXO J013329.0+304216 


1730 


bremss 


0.34+«,f 


^•^•^-0.08 






786 


bremss 


S7+0-44 


'^•^^-0.09 






2023 


bremss 


ti5+0-O4 
^J-OO-O.Ol 


'-'•i"-o.o 












kT 


z 




1730 


apec 


< 0.31 


0.36+!5f 


11+019 
^•11-0.07 




786 


apec 


< 0.16 


4+0.14 


'^•^^-0.12 




2023 


apec 


0.25l0;i5 


0.28+0.05 












kT 


z 


CXO J013329.3+304508 


1730 


apec 


14+°-^ 


^•'-1.4 


< 4.5 




786 


apec 


< 0.13 


o r, + 16.8 


< 1.9 




2023 


apec 


< 0.09 


1 7 o + oo 
l'-^-10.2 


2.0" 










kT 




CXO J013329.4+304912 


786 


bremss 


'^■^^-0.34 


^•l>'-0.07 






2023 


bremss 


,-+0.08 
-0.02 


n 9+0.08 
^•^-0.01 












kT 


Z 




786 


apec 


< 0.31 


0.3ll°f 


fi Q+2.2 
"■•^-0.26 




2023 


apec 


< 0.6 


26+°-*^ 
U.^O_0.13 


f, nr;+0.23 



kT Ei/2 0-1/2 

CXO J013331.1+303333 1730 bremss 0.17t'^ °l 0.291°;°^ 0.89lg < 0.045 

(M33 X-14) +2*gauss 1.35lg;°2 < 0-035 
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TABLE Bl — Continued 



Source 


Obsid 


XSPEC model 


Nh 
[1022 cm-2] 


parameters 
kT, TO, Tin, E, a: [keV], 
Z: [solar abund.] 










Tin 




CXO J013334. 1+303210 
(M33 X-7) 


1730 


diskbb 


< 0.05 


1 n+O-05 

^•^-0.03 

kT 




CXO J013335.9+303627 


1730 
1730 


bremss 
apec 


< 0.38 


si+0" 

kT 
0.40+°?| 

— U. iz 


z 

0.12+°-5B 
— u.uy 


CXO J013336.3+303742 


1730 
786 


pow 
pow 


62+°-'^ 
'-'•DZ_Q 33 


r 

9 1+0.9 

^•1-0.8 
1 9+0.85 
-^•^-0.76 

r 




CXO J013337.4+304718 


1730 
786 
2023 


pow 
pow 
pow 


< 0.33 
n 22+°-" 

< 0.49 


1 Q+0.65 
-'-■^-0.55 
n+0-55 
^•^-0.45 

n s+0-7 

0-°-0.66 

kT 


z 


CXO J013341. 8+303848 


1730 
786 
2023 


apec 
apec 
apec 


< 0.07 

< 0.05 

< 0.07 


0.70±0.15 

(, -„+0.09 

(, -9+0.17 
^•'^-0.15 

kT 


0.17+°}^ 
19+0.07 
"■-'-''-0.06 

< 0.04 


CXO J013343.4+304630 


1730 

786 
2023 


bbody 

bbody 
bbody 


0.79* 

0.79l°-t2 

0.83l0:g2 


0.07* 

r, r,7 + 0.035 

0-OS_o.055 

K i 


K i Zj 


CXO J013346.2+303807 


2023 


bbody+apec 


0.06* 


0.02±0.01 


3.lt.ll 0.3* 


CXO J013346.5+303748 


1730 
786 


pow 
pow 


< 0.14 


r 

1.6±0.35 
1 1+0.35 

^■^-0.26 

r 




CXO J013350.5+303821 


1730 
786 


pow 
pow 


< 0.39 

< 0.16 


1 fi+O-S 
l-D-0.7 
1 r;+0.4 
J-'->_0.4 

kT 




CXO J013354.8+303309 
(M33 X-13) 


1730 
786 


bremss 
bremss 


o.3oi?,:?^ 


99+0.07 
"J-^^-0.07 
-,0+0.04 
0-1°-0.05 




CXO J013356.8+303729 


1730 
786 


pow 
pow 


0.15* 


r 

1.5* 

1 5+°-* 
1-^-0.36 

r 




CXO J013358.8+305004 
(M33 X-9a) 


786 
2023 


pow 
pow 


< 0.46 
54+°-22 


-1 7+0.2 
-0.8 

1 S+O-5 

kT 


z r 


CXO J013401.5+303136 


786 


apec+pow 


n Q4+0-5 
"•^^-0.54 


6.3" 
kT 


0.3* 3.4l5;2 

r 


CXO J013402.0+303004 


786 


bbody+pow 


< 0.27 


kT 


1 0+0.86 

^•**-0.35 


CXO J013409.9+303219 


786 


bbody 


70+0.34 


(, „o~+0.007 
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TABLE Bl — Continued 



Source 


Obsid 


XSPEC model 


Nh 
[1022 cm-2] 


parameters 
kT, TO, Tin, E, a: [keV], 
Z: [solar abund.] 


CXO J013410.5+303946 


1730 


diskbb 


0.06^ 


1 2^^^- 
^•^-0.3 






786 
2023 


pow 
diskbb 


0.06* 
0.06* 


r 

1.7±0.1 

Tin 
1 2it0 2 

kT 




CXO J013410.6+304223 


1730 


bremss 


U. 'y_o.43 


11+00^ 






2023 


bremss 


< 0.19 


n Q+0.07 

Tin 




CXO J013421.2+304932 


2023 


diskbb 


1 17+0.55 


1 7+1-'' 
-^•'-0.6 


^break ^2 


CXO J013425. 4+302821 


786 


bknpow 


< 0.25 




> 1.2 o.ieliio 


CXO J013425.5+305514 


786 
2023 


pow 
pow 


17+0.03 
"•^'-0.04 


r 

r, 1+0.25 

^•-^-0.15 
^■'-'-0.04 

i 




CXO J013427.0+304314 


786 


pow 


< 0.27 


9+1-6 
"•''-0.3 






2023 


pow 


< 0.17 


^-1-0.6 
^•'-0.3 

kT 




CXO J013432.0+303455 


1730 


bremss 


27+0.8 


r, y+11.3 

^- ' -0.9 

kT 


r 


CXO J013433.7+304701 


786 
2023 


bbody+pow 
bbody+pow 


< 0.9 


O.liO.03 

o.o9+;;-!S 

^ — 0.04 

kT 
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Fig. B1.— 



M33 X-ray variability 

X-ray spectra and contour plots for best fit model. Models and parameters are shown in Table Bl. 
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